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FOREWORD 


This  report  documents  an  investigation  of  the  electrical, 
thermal,  and  thermostructural  performance  of  graphite  heating 
elements  used  in  the  Naval  Surface  Warfare  Center  (NSWC) 
Hypervelocity  Wind  Tunnel  No.  9.  This  study  is  part  of  a  larger 
and  ongoing  effort  by  the  Aerodynamic  Facilities  Branch 
(Code  K23)  to  reduce  or  eliminate  the  problem  of  frequent  and 
costly  failures  of  the  graphite  heaters.  The  analysis  herein 
focuses  on  the  nitrogen  preheating  period  of  the  tunnel  run  cycle 
during  which  most  heater  failures  are  known  to  occur. 

A  general  finite  element-based  method  for  investigating  the 
performance  of  electrically  powered  heating  elements  is 
presented.  The  technique  is  applied  to  the  Tunnel  No.  9  graphite 
heater,  and  it  is  shown  to  correctly  predict  brittle  fracture  at 
fillet  locations  in  the  heater  body  where,  in  fact,  most  heater 
fractures  occur.  Two  methods  for  reducing  fillet  stresses  are 
proposed  which,  if  implemented,  could  substantially  increase  the 
useful  service  life  of  the  graphite  heaters. 

The  author  wishes  to  thank  Dr.  Robert  Edwards  and  the  K22 
Structural  Mechanics  Group  members  for  their  guidance  concerning 
this  study.  Also,  thanks  are  extended  to  Margie  Fung,  Aero- 
Thermodynamics  Group  (Code  K22) ,  for  providing  the  procedure  for 
calculating  free-convection  film  coefficients,  to  Mr.  Ray 
Trohanowsky  for  developing  the  POWEF  -tran  subroutine.  And, 
finally,  thanks  to  the  author's  brai.ch  head  Dr.  J.  Michael 
Etheridge  for  the  many  useful  suggest  -  s  for  improving  this 
technical  report. 


Accesion  For 

NTIS  CRA&I 
OTIC  TAB 

Unannounced 

Justification 


By _ 

Distribution/ 


J 

a 

□ 


Availability  Codes 


Dist 


M 


Avail  unrj/or 


Approved  by: 


L.  SCHMIDT,  Head 
Strategic  Systems  Department 


iii/iv 


NSWC  TR  88-146 


CONTENTS 

Chapter  Page 

1  INTRODUCTION  .  1 

BACKGROUND  .  1 

TUNNEL-9  OPERATION  .  1 

HEATER  ELEMENT  BREAKAGE  PROBLEM . 2 

2  GENERAL  HEATER  ELEMENT  PERFORMANCE  ANALYSIS 

PROCEDURE  .  5 

ANALYSIS  PROCEDURE  .  5 

ELECTRIC  RESISTANCE  HEATING — STEP  1 .  6 

GRAPHITE  ELEMENT  HEATING— STEP  2 . 7 

GRAPHITE  ELEMENT  THERMAL  STRESSES— STEP  3.  .  .  8 

PRACTICUM . 8 

3  PERFORMANCE  ANALYSIS  RESULTS  FOR  TUNNEL-9 

GRAPHITE  HEATING  ELEMENTS  .  11 

ANALYSIS  RESULTS  .  11 

ELECTRICAL  ANALYSIS  RESULTS— STEP  1 .  11 

THERMAL  ANALYSIS  RESULTS— STEP  2 . 12 

THERMOSTRUCTURAL  ANALYSIS  RESULTS — STEP  3 .  .  .  12 

GRAPHITE  HEATER  ELEMENT  FILLET  STRESSES  .  15 

FILLET  STRESS  BREAKDOWN . 15 

STRESS  REDUCTION  METHODS  .....  .  16 

ASSUMPTIONS . 17 

4  CONCLUSIONS  .  18 

REFERENCES . 63 

DISTRIBUTION  .  (1) 


&PP£ndi.X  Page 

A  ANALYSIS  SUPPORT  PROGRAMS . A-l 

B  ABAQUS  INPUT  FILES  .  B-l 

C  FREE  CONVECTION  HEAT  TRANSFER  FILM  COEFFICIENT 

CALCULATION  FOR  MACH-14  END  HEAT  CONDITIONS.  .  .  C-l 

D  TENSILE  STRENGTH  AND  STRAIN-TO-FAILURE 

EXPERIMENTAL  DATA  FOR  HEATER  GRAPHITE . D-l 


V 


NSWC  TR  88-146 


ILLUSTRATIONS 

Figure  Page 

1  SCHEMATIC  OF  TUNNEL  NO.  9  FACILITY . 19 

2  VERTICAL  HEATER  AND  DIAPHRAGM  SECTION  .  19 

3  HEATER  GAS  CONDITIONS  FOR  A  TYPICAL  MACH- 14 

WIND  TUNNEL  RUN . 20 

4  COMPLETELY  FAILED  HEATER — TOP  FILLET  LOCATION  ...  21 

5  COMPLETELY  FAILED  HEATER— BOTTOM  FILLET  LOCATION.  .  21 

6  TOP  AND  BOTTOM  FILLET  CRACK  LOCATIONS  IN  THE 

GRAPHITE  HEATER . 22 

7  PARTIAL  TOP  FILLET  CRACK . 23 

8  PARTIAL  BOTTOM  FILLET  CRACKS . 23 

9  HEATER  GRAPHITE  ELECTRICAL  RESISTIVITY 

VERSUS  TEMPERATURE  .  24 

10  J  HEATER  DRAWING . 25 

11  K  HEATER  DRAWING . 26 

12  MACH-10  HEATER  BASE  DRAWING . 27 

13  J  HEATER  FINITE  ELEMENT  MODEL . 28 

14  K  HEATER  FINITE  ELEMENT  MODEL . 28 

15  FINITE  ELEMENT  MODEL  USES  HEATER  SYMMETRY  TO 

REDUCE  MODEL  SIZE . 29 

16  MACH-10  HEATER  BASE  FINITE  ELEMENT  MODEL . 29 

17  VOLTAGE  POTENTIALS  IN  THE  J  HEATER . 3  0 

18  CURRENT  DENSITY  CONTOURS  IN  THE  J  HEATER . 31 

19  START  HEAT  TEMPERATURES  IN  J  HEATER . 3  3 

20  END  HEAT  TEMPERATURES  IN  THE  J  HEATER . 33 

21  END  HEAT  TEMPERATURES  IN  THE  J  HEATER . 35 

22  START  HEAT  TEMPERATURES  IN  THE  K  HEATER . 35 

23  START  HEAT  TEMPERATURES  IN  THE  MACH-10  BASE  ....  37 

24  AXIAL  STRESS  IN  MACH-10  BASE . 37 

25  HEATER  ELEMENT/BASE  CLAMPING  ARRANGEMENT . 39 

26  DEFORMED  SHAPE  PLOTS —FLOATING  BASE . 40 

27  DEFORMED  SHAPE  PLOTS— CLAMPED  BASE . 41 

28  AXIAL  STRESS  IN  J  HEATER . 43 

29  AXIAL  STRESS  IN  K  HEATER . 4  3 

30  MAXIMUM  PRINCIPAL  STRESS — J  HEATER,  START  HEAT.  .  .  45 

31  MAXIMUM  PRINCIPAL  STRESS — J  HEATER,  END  HEAT.  ...  45 

32  MAXIMUM  PRINCIPAL  STRESS— J  HEATER,  FREE  BASE  ...  47 

33  MAXIMUM  PRINCIPAL  STRESS— K  HEATER,  START  HEAT.  .  .  47 

34  MAXIMUM  PRINCIPAL  STRESS— MACH  10  BASE . 49 

35  TEMPERATURES  IN  SQUARE  AND  ROUNDED  CORNER 

HEATER  LEG . 49 

vi 


NSWC  TR  88-146 


ILLUSTRATIONS  (Cont.) 

Figure  Paqe 

36  TOTAL  PRINCIPAL  STRESS  AND  TEMPERATURE  PROFILES 

EXTRAPOLATED  TO  HEATER  SURFACE. 

TOP  FILLET  LOCATION . 51 

37  TOTAL  PRINCIPAL  STRESS  AND  TEMPERATURE  PROFILES 

EXTRAPOLATED  TO  HEATER  SURFACE. 

BOTTOM  FILLET  LOCATION  .  52 

38  TOTAL  PRINCIPAL  ELASTIC  STRAIN  PROFILE  EXTRAPOLATED 

TO  HEATER  SURFACE.  TOP  FILLET  LOCATION . 53 

39  TOTAL  PRINCIPAL  ELASTIC  STRAIN  PROFILE  EXTRAPOLATED 

TO  HEATER  SURFACE.  BOTTOM  FILLET  LOCATION  ....  53 

40  GRAPHITE  HEATER  SERVICE  LIFE  IN  TUNNEL  RUNS  VS. 

CUMULATIVE  FREQUENCY . 54 

41  BREAKDOWN  OF  THERMAL  STRESS  AT  BOTTOM  FILLET.  ...  55 

42  FINITE  ELEMENT  MODEL  WITH  SLOT  EDGES  ROUNDED  TO 

REDUCE  FILLET  STRESSES . 56 

43  COOLING  SLOTS  FORMED  IN  HEATER  LEGS  BY  SECTIONING 

AT  AA  AND  SEPARATING  HALVES . 56 

44  CONCEPT  FOR  HEATER  ELEMENT  WITH  COOLING  SLOTS  ...  57 


TABLES 

Table  Page 

1  HEATER  ELEMENT  SERVICE  LIFE  HISTORY  FOR  1985 

AND  1986 .  58 

2  SUMMARY  OF  HEATER  ANALYSIS  STEPS . 59 

3  HEAT  TRANSFER  INPUT  DATA . 60 

4  PEAK  THERMAL  STRESS -STRAIN  STATES  IN  GRAPHITE 

HEATER . 61 

5  PROBABILITY  OF  FAILURE  FOR  GRAPHITE  HEATERS  ....  61 


vii/viii 


NSWC  TR  88-146 


CHAPTER  1 
INTRODUCTION 


The  Naval  Surface  Warfare  Center  (NSWC)  Hypervelocity  Wind 
Tunnel  No.  9  uses  an  electrically  powered  graphite  heating 
element  to  preheat  nitrogen  working  gas  to  temperatures  up  to 
3100°F.  This  prevents  condensation  of  the  gas  which  cools  as  it 
expands  through  a  nozzle  to  hypersonic  speeds  of  Mach  10  or 
Mach  14.  The  graphite  heater  elements  have  a  limited  and  erratic 
service  life — a  newly  installed  heater  element  may  fail  after 
only  1  or  as  many  as  100  or  more  tunnel  run  cycles.  Element 
replacement  is  costly  (approximately  $10,000  including 
installation) ,  and  a  half  day  or  more  may  be  required  to  replace 
a  broken  heater  element  which  adversely  impacts  the  tunnel  test 
schedule.  This  report  documents  a  general  method  for 
investigating  the  electrical,  thermal,  and  thermostructural 
performance  of  electrically  powered  heating  elements.  It  also 
presents  the  results  of  an  investigation  of  the  heating 
characteristics  of  the  Tunnel-9  graphite  heater  elements.  This 
study  was  initiated  to  determine  if  thermal  stresses  induced  in 
the  heater  elements  during  the  nitrogen  preheating  period  could 
be  the  cause  of  the  heater  element  failures. 


BACKGROUND 
T^nn^lr?  .Operation 

The  NSWC  Tunnel-9  wind  tunnel  is  a  "blow-down"  type  tunnel 
which  runs  at  Mach  10  or  14.  High  pressure  vertical  heater  and 
driver  vessels  (shown  in  Figure  l)  are  precharged  with  nitrogen 
gas  prior  to  a  wind  tunnel  run  while  a  vacuum  sphere  downstream 
of  the  test  section  is  evacuated.  The  fixed  volume  of  gas  in  the 
heater  vessel  is  preheated  to  the  extreme  temperature  required. 

A  tunnel  run  is  initiated  by  bursting  a  set  of  metal  diaphragms 
located  just  upstream  of  a  40-foot  long  nozzle  which  allows  the 
high-pressure,  high-temperature  nitrogen  (typically  20,000  psi  at 
3100°F  for  Mach  14)  to  accelerate  through  the  nozzle  section  and 
reach  speeds  of  Mach  10  or  14  in  the  5-foot  diameter  test 
section.  Relatively  cool  driver  vessel  gas  pushes  the  hot  gas 
out  of  the  heater  vessel  under  constant  pressure  to  maintain 
constant  flow  conditions  in  the  test  section  during  the  tunnel 
run.  The  tunnel  run  continues  until  all  of  the  hot  gas  has  been 
expended,  a  typical  run  lasting  approximately  from  .25  second  to 
as  long  10  seconds. 
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The  heater  vessel  (shown  in  Figure  2)  is  lined  with  a 
refractory  insulation  liner  which  safely  contains  the  nitrogen 
gas  during  the  preheating  process.  The  liner  encloses  the 
graphite  heater  element  which  heats  the  nitrogen  precharge  to  the 
specified  final  temperature  of  1500°F  (for  Mach  10)  or  3100°F 
(for  Mach  14).  This  heating  process  is  depicted  in  Figure  3. 
During  the  heating  period,  which  lasts  about  15  minutes,  a  60  Hz, 
single  phase  ac  current  of  5500-6500  rms  amps  is  passed  through 
the  hairpin-shaped  graphite  heater  element  which  reaches  a 
temperature  of  4000-5000°F  as  a  result  of  resistance  heating. 

The  nitrogen  gas,  which  is  essentially  transparent  to  radiation, 
is  heated  primarily  by  convection  and  conduction  off  the  hot 
graphite  element  and  hot  liner  walls.  Since  the  volume  of 
nitrogen  in  the  heater  vessel  is  fixed,  both  its  temperature  and 
pressure  rise  during  the  constant  volume  heating  process  until 
the  final  desired  pressure  and  temperature  are  reached.  At  this 
point  the  tunnel  flow  is  started  by  rupturing  the  metal  sealing 
diaphragms . 


Heater  Element  Breakage  Problem 

Service  life  histories  for  heater  elements  used  during  FY 
1985  and  1986  are  shown  in  Table  1.  A  heater  element  is 
typically  operated  until  it  fails  " catastrophically, "  usually  by 
fracturing  completely  through  one  of  the  heating  legs  as  shown 
in  Figures  4  and  5.  The  timing  and  mode  of  failure  are  fairly 
consistent.  The  element  usually  fails  completely  at  some  point 
during  the  nitrogen  heating  period.  The  fracture  is  always  a 
brittle-type  fracture  typically  located  at  either  a  top  fillet  or 
bottom  fillet  (as  shown  in  Figure  6)  with  the  crack  surfaces 
oriented  approximately  perpendicular  to  the  longitudinal  axis  of 
the  heater  element. 

Interestingly,  it  is  now  fairly  well  documented  that  a 
heater  may  develop  partial  cracks  at  the  top  or  bottom  fillets 
(as  shown  in  Figures  7  and  8)  which,  at  least  for  a  time,  will 
not  adversely  affect  the  heater  operation.  However,  eventually 
the  element  fails  completely  and  must  be  replaced.  Although 
complete  element  failure  occurs  during  the  heating  period,  it  is 
not  clear  if  a  partial  crack  is  initiated  at  this  or  some  other 
time. 


Other  work  has  been  done  to  try  to  improve  heater  element 
life.  Results  of  an  investigation  of  thermal  stresses,  induced 
by  quenching  of  the  heater  as  the  cool  driver  gas  fills  the 
heater  vessel  during  a  tunnel  run,  did  not  correlate  well  with 
Tunnel-9  experience.  In  another  study,  graphite  material  samples 
taken  from  both  long-lived  and  short-lived  failed  heater  elements 
were  tested  in  an  attempt  to  correlate  service  life  with  certain, 
easily  obtained,  non-destructive  test  measurements  which  could 
then  be  used  as  a  quality  control  criterion  for  screening  out 
unacceptable  graphite  material  billets;  however,  no  definite 
correlations  were  found.1 
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The  present  study  is  concerned  with  determining  if  thermal 
stresses  induced  in  the  heater  elements  during  the  nitrogen 
preheating  period  could  be  the  cause  of  heater  element  failures. 
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CHAPTER  2 

GENERAL  HEATER  ELEMENT  PERFORMANCE  ANALYSIS  PROCEDURE 


ANALYSIS  PROCEDURE 

The  heater  analysis  procedure  can  be  broken  down  into  three 
basic  steps: 

Step  1:  Solve  the  Electrical  Problem — Get  volumetric 
resistance  heating  in  the  graphite  heater 
generated  by  a  (5500  amp)  current  flow  through  the 
heater . 

Step  2:  Solve  the  Thermal  Problem — Use  internal  heating 

from  Step  1  to  ’’heat  up"  the  graphite  element  and 
obtain  graphite  temperatures. 

Step  3 :  Solve  the  Structural  Problem 

o  Use  temperatures  from  Step  2  to  obtain  thermal 
stress  and  strain  in  heater  due  to  thermal 
expansion.  Determine  maximum  principal  stress 
and  strain  in  the  heater. 

o  Subtract  the  hydrostatic  stress  or 

strain  (generated  by  the  high-pressure  nitrogen 
gas  surrounding  the  heater)  from  the  peak 
principal  stress  or  strain  to  obtain  the  net 
maximum  principal  stress  or  strain. 

o  Compare  net  maximum  principal  stress  or  strain 
with  graphite  strength  to  determine  the 
probability  of  failure. 


The  entire  three-step  procedure  above  was  implemented  using 
the  ABAQUS2  Finite  Element  Analysis  software  (version  4-5)  in 
conjunction  with  the  PATRAN3  pre-  and  post-processing  code.  The 
implementation  of  each  of  the  three  steps  is  described  in  more 
detail  below. 
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Electric  Resistance  Heating— Step  1 

Implementation  of  Step  1,  that  is,  obtaining  the  volumetric 
resistance  heating  (also  known  as  joule  heating)  due  to  the 
electric  current  passing  through  the  graphite  heater  element, 
required  a  novel  application  of  the  well  documented 
thermal/electric  (T/E)  analogy  principal.4  The  T/E  analogy  has 
been  widely  used  to  solve  analogous  engineering  problems  in  the 
fields  of  heat  conduction,  fluid  dynamics,  and  solid  mechanics. 
Phenomena  which  are  governed  by  the  Laplace  equation,  such  as 
both  electrical  and  heat  conduction,  yield  mathematically 
identical  solutions;  the  solution  of  a  particular  electrical 
conduction  problem  is  also  applicable  to  the  corresponding  or 
analogous  heat  conduction  problem  and  vice  versa.  Analogous 
mathematical  terms  for  the  T/E  analogy  are  given  below: 


THERMAL 


Heat,  Q 
Temperature,  T 
Heat  Flux  Vector,  q 
Thermal  Conductivity,  k 


ELECTRICAL 


Charge ,  Q 
Voltage,  V 

Current  Density  Vector,  j 
Electrical  Conductivity, s 


The  typical  approach  to  applying  the  T/E  analogy  has  been  to 
obtain  a  solution  to  a  particular  heat  conduction  problem  using 
an  analogous  electrical  conduction  model  (e.g.,  electrolytic 
tanks  or  conductive  sheet  models) .  In  the  present  analysis, 
however,  the  opposite  approach  has  been  taken;  that  is,  the 
electrical  conduction  in  the  heater  is  obtained  by  solving  the 
analogous  heat  conduction  problem  using  a  geometrically  similar 
finite  element  heat  conduction  model. 


The  ABAQUS  thermal  conduction  capability  can  be  utilized  to 
easily  solve  for  current  flow  in  an  arbitrarily  shaped  conductor 
subjected  to  some  prescribed  voltage.  Since  ABAQUS  does  not 
"know"  that  it  is  solving  an  analogous  electrical  conduction 
problem,  it  will  state  the  results  as  "temperatures "  and  "heat- 
flux,"  etc.  However  if  consistent  analogous  electrical  terms  are 
input,  the  "temperature"  and  "heat-flux"  results  from  ABAQUS  will 
actually  be  voltage  and  current-density,  respectively.  For  this 
analysis,  the  units  used  for  the  input  quantities  were  inches, 
volts  and,  for  electrical  conductivity,  1/ohm-in  which  produce 
units  of  volts  and  amps/sq  in  for  the  voltage  and  current- 
density  output  quantities, respectively.  ABAQUS  will  calculate 
the  components  and  also  the  magnitude  of  the  current-density 
vector,  j,  at  each  gauss  point  in  the  finite  element  model.  The 
volumetric  resistance  heating  at  each  gauss  point  can  then  be 
obtained  by  Equation  (1) : 


q  -  |j|2P  U) 

where:  q  «  volumetric  heat  flux  ( watts/ cu  in) 

| j |  -  current  density  magnitude  (amps/sq  in) 
p  ■  electrical  resistivity  (ohm-in) 
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This  last  step  of  calculating  the  heat  flux  at  each  gauss 
point  is  done  by  running  a  short  Fortran  program  I2RFIL  (listed 
in  Appendix  A) .  This  program  simply  reads  the  current  density 
values  from  the  ABAQUS  results  (FIL)  file  and  generates  a  file  of 
volumetric  heat  flux  at  each  gauss  point  using  Equation  (1) . 

This  file  of  heat  flux  values  is  then  used  to  heat-up  the 
graphite  element  (Step  2  in  the  analysis  procedure)  to  obtain  the 
temperature  in  the  heater. 

Limitations  of  T/E  Method.  There  are  two  principal 
limitations  in  the  T/E  method  described  above.  First,  the 
electrical  conductivity  must  remain  constant  with  temperature. 
Although  ABAQUS  will  handle  temperature  dependent  thermal 
conductivity,  this  capability  cannot  presently  be  utilized  for 
temperature  dependent  electrical  conductivity  since  the  actual 
temperatures  are  not  known  in  Step  1.  They  are  computed  later  on 
during  Step  2  in  the  analysis  procedure.  As  Figure  9  shows,  the 
electrical  resistivity  (the  resistivity  is  simply  the  inverse  of 
conductivity)  of  the  graphite  material  is  strongly  temperature 
dependent  up  to  about  1400°F,  then  remains  fairly  constant  at 
temperatures  above  this.  It  was  found  that  the  heater  element 
legs  run  at  temperatures  above  this,  except  for  the  initial  heat¬ 
up  period,  which  lasts  about  a  minute,  so  that  for  most  of  the 
15-20  minute  nitrogen  heating  period  the  electrical  conductivity 
can  be  taken  to  be  uniform  and  constant. 

This  method  is  also  in  a  strict  sense  limited  to  dc 
electrical  power.  Recall  that  the  Tunnel-9  heater  elements  are 
powered  by  a  60  Hz,  ac  power  source.  An  ac  current  produces  a 
so-called  "skin  effect"  which  causes  a  concentration  of  the 
current  density  at  the  outer  surface  of  a  conductor  rather  than 
being  uniformly  distributed  through  the  cross  section  as  in  dc 
current  flow.  This  effect  is  proportional  to  the  ac  frequency 
and  size  of  the  conductor  and  is  inversely  proportional  to 
resistivity.  It  was  determined  that,  in  regard  to  the  volumetric 
heating,  this  effect  would  be  secondary  for  the  graphite  heater 
elements  so  that  the  dc  current  flow  analogy  should  adequately 
model  the  heating  characteristics. 


graphite  .El K3atenq-“?tep  z 

Since  both  Steps  1  and  2  in  this  analysis  procedure  use  a 
heat  condu c*  >n  finite  element  model,  the  basic  model  used  in 
Step  1  can  .  ^  used  in  Step  2  with  changes  in  the  boundary 
conditions  l  reflect  the  new  nature  of  the  problem.  ABAQUS  is 
capable  on  solving  very  general  transient  and  steady-state  heat 
transfer  pro*''  ts  involving  internal  heat  generation,  radiation, 
and  conv  ctio.,  with  temperature  dependent  properties. 

Recall  from  Figure  3  that  once  power  to  the  heater  is  turned 
on  at  the  start  of  the  heating  period,  the  graphite  element  heats 
up  in  about  1  minute.  During  this  time  the  temperature  of  the 


7 


NSWC  TR  88-146 


nitrogen  gas  in  the  heater  vessel  does  not  increase  appreciably. 
For  analysis  purposes  this  initial  heating  transient  must  be 
viewed  as  a  nonsteady-state  phenomenon.  After  about  1  minute, 
the  heater  element  temperature  reaches  a  nearly  steady-state 
level  and,  thereafter,  both  the  heater  and  nitrogen  gas 
temperatures  rise  slowly  for  the  remaining  15-20  minutes  of 
heating  time.  For  analysis  purposes  this  relatively  slow  heating 
process  can  be  treated  as  a  quasi-steady  state  process  so  that 
the  heater  temperature  can  be  assumed  to  be  a  function  solely  of 
the  surrounding  gas/liner-wall  temperature.  The  problem  of 
obtaining  the  temperature  distribution  in  the  heater  at  any  point 
during  the  quasi-steady  heating  reduces  to  simply  fixing  the 
gas/wall  temperature  and  solving  for  the  steady-state  heater 
temperature  distribution.  Since  performing  a  transient  analysis 
is  in  general  much  more  costly,  it  was  decided  to  limit  the 
present  study  to  investigating  the  heater  temperatures  and 
stresses  at  selected  points  in  time  within  this  quasi-steady 
heating  period.  The  two  points  selected  for  study  were  at  the 
start  and  end  of  the  quasi-steady  heating  period  and  which  are 
designated  as  "START  HEAT"  and  "END  HEAT"  in  Figure  3. 
Additionally,  the  gas  and  liner-wall  temperatures  were  assumed  to 
be  identical  permitting  the  use  of  one  sink  temperature  for  both 
radiation  and  convection  boundary  conditions. 

Obtaining  the  temperature  distribution  in  the  graphite 
heater  at  the  START  or  END  HEAT  is  a  matter  of  specifying  the 
proper  sink  temperature  in  ABAQUS,  as  well  as  specifying 
radiation  and  convection  heat  transfer  coefficients,  and  then 
allowing  ABAQUS  to  solve  for  the  steady-state  temperatures  in  the 
heater  due  to  the  resistance  heating  provided  by  Step  1.  Since 
the  resistance  heating  in  the  heater  is  non-uniform,  the  ABAQUS 
subroutine  option  DFLUX  is  used  with  the  subroutine  given  in 
Appendix  A  to  read  the  heating  flux  values  from  the  file 
generated  in  Step  1.  ABAQUS  calculates  the  temperatures  in  the 
heater  which  are  then  used  to  evaluate  thermal  performance  of  the 
heater  and  to  obtain  thermal  stresses  in  the  heater. 


Graphite  Element  Thermal  Stresses — Step  3 

Once  again  the  same  basic  finite  element  model  from  Steps  1 
and  2  can  be  used  in  Step  3  if  the  element  type  designator  is 
changed  to  the  corresponding  solid  structural  element.  ABAQUS 
will  then  use  the  temperatures  obtained  from  Step  2  to  compute 
the  thermal  stresses  and  strains  in  the  graphite  heater  element. 
A  suitable  failure  criteria  can  then  be  applied  to  the  stress  or 
strain  results  to  determine  probability  of  failure  or  assess 
margins  of  safety  of  the  heater  design. 


PRACTICUM 

In  practice,  implementation  of  this  three-step  analysis 
procedure  involves  additional  work,  primarily  in  Step  1,  to 
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ensure  that  correct  results  are  obtained.  Table  2  summarizes  the 
semi-automated  steps  for  the  entire  procedure  as  implemented  on  a 
VAX  computer  system.  In  Step  1  the  user  specifies  a  heater 
voltage  and  ABAQUS  solves  for  current  flow.  The  current  checking 
program  CURRENT  (Appendix  A)  will  verify  that  the  net  current 
flow  through  the  heater  is  correct.  If  it  is  not,  then  the 
applied  voltage  may  be  adjusted.  Once  the  correct  net  current  is 
obtained,  the  program  I2RFIL  (Appendix  A)  is  run  which  creates 
the  resistance  heating  file  (FOR032.DAT  file)  from  the  ABAQUS 
results  file.  Once  the  heating  file  is  generated,  it  should  be 
checked  by  running  the  POWER  program  (Appendix  A)  which 
calculates  the  total  heater  power  that  will  be  generated  by  the 
heating  file.  The  total  power  calculated  by  POWER  should  be 
within  a  few  percent  of  the  specified  heater  voltage  multiplied 
by  the  total  heater  current  computed  by  CURRENT  ( i . e . ,  the  total 
resistance-heating  power  should  equal  applied  heater  volts  times 
total  sumps)  .  Some  discrepancy  can  be  expected  since  the 
algorithm  used  in  POWER  to  compute  finite  element  volumes  is  not 
exact  for  quadratic  elements  with  curved  edges. 
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CHAPTER  3 

PERFORMANCE  ANALYSIS  RESULTS  FOR  TUNNEL-9  GRAPHITE 

HEATING  ELEMENTS 


ANALYSIS  RESULTS 

The  two  most  current  graphite  heater  element  configurations, 
hereafter  referred  to  as  the  J-type  and  K-type  heaters,  are  shown 
in  Figure  10  and  11.  The  drawing  for  the  Mach  10  heater  base  is 
shown  in  Figure  12.  The  K  heater  has  superceded  the  J  heater, 
the  differences  being  in  the  flared  top  and  the  side  slots.  In 
the  K  design,  the  slight  flare-out  of  the  outer  diameter  at  the 
top  end  was  lengthened,  and  the  slot  was  also  lengthened  in  order 
to  "beef-up"  the  top  fillet  region  where  the  earlier  J-type 
heaters  tended  to  break.  This  change  appears  to  have  reduced 
top  fillet  breaks  but,  unfortunately,  the  K  heaters  continued  to 
fracture  at  the  bottom  fillet.  Both  the  J  and  K  heaters  were 
analyzed  with  the  basic  finite  element  models  shown  in  Figure  13 
and  14.  Both  are  l/4-symmetry  models  utilizing  the  two  planes  of 
symmetry  in  the  heater  to  reduce  model  size  as  shown  in  Figure 
15.  The  J  model  is  a  full-length  model,  whereas  only  the  top 
portion  of  the  K  design  was  analyzed.  Extensive  studies  showed 
that  such  top  models  would  give  the  same  stress  results  as  the 
full  model  when  out-of-plane  motion  of  the  cross  section  at  the 
cutoff  is  prevented.  Quadratic  finite  elements  were  used  in  both 
the  thermal  models  and  the  stress  models  (20-node  hexes  and  15- 
node  wedges,  ABAQUS  DC3D20-  and  C3D20-type  elements, 
respectively) .  The  Mach  10  graphite  heater  base,  which  holds  the 
heater  element,  was  also  analyzed  separately  using  the  1/4- 
symmetry  model  shown  in  Figure  16.  Only  the  bottom  portion  of 
the  heater  element,  up  to  the  bottom  fillet  region,  is  included 
in  the  base  model.  The  results  of  the  base  analysis  will  not  be 
discussed  extensively  in  this  report— these  results  are  presented 
for  reference. 


Electrical  Analysis  Results— Step  1 

Results  for  voltage  and  current  density  distributions  in  the 
J-type  heater  are  shown  in  Figure  17  and  18,  respectively.  The 
current  density  plot  can  be  viewed  as  a  topographical  map  with 
peaks  showing  regions  of  high  current  densities,  such  as  at  both 
the  top  and  bottom  fillets,  and  valleys  indicating  areas  of 
relatively  low  current  densities.  The  ABAQUS  material  and 
boundary  conditions  are  given  in  Appendix  B  for  the  J  heater 
electrical  model. 
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Thermal  Analysis.  Results— Step_  2 

Thermal  Model  Setup.  In  all  cases  convection  heat  transfer 
was  allowed  on  the  inner  diameter  (I.D.)  and  outer  diameter 
(0.0.)  "wetted"  surfaces  of  the  heater  element  as  well  as  on  the 
rounded  end  of  the  crown  and  edges  of  the  side  slots.  The  O.D. 
and  I.D.  of  the  threaded-base  portion  of  the  J  heater  was  fixed 
at  300°F.  Radiation  heat  transfer  was  allowed  on  the  wetted 
O.D. ,  slot  edges,  and  the  rounded  end  of  the  crown.  No  radiation 
was  allowed  on  the  wetted  I.D.  since  the  inside  surfaces  of  the 
heater  element  legs  face  each  other  so  that  the  net  radiation  to 
each  surface  is  taken  to  be  zero  as  a  first  approximation.  Table 
3  lists  the  film  coefficients  and  radiation  boundary  conditions 
used.  ABAQUS  input  data  for  the  J  thermal  model  are  given  in 
Appendix  B,  and  Appendix  C  shows  the  method  used  to  compute  the 
free  convection  film  coefficient. 

Thermal  Results.  Figures  19-23  show  temperature  results  for 
the  three  models  analyzed.  Figures  19  and  20  show  temperatures 
in  the  J  heater  for  a  standard  613  kW  power  case  (5481.8  A  @  112 
V)  for  both  START  HEAT  and  END  HEAT  periods.  The  peak 
temperatures  in  the  graphite  are  3620°F  at  START  HEAT  and  4858°F 
at  END  HEAT,  and  they  occur  at  the  "hotspot"  apparent  on  the  I.D. 
surface  near  the  top  fillet.  Cross-sectional  temperature 
profiles  at  the  top  and  bottom  fillets  for  the  END  HEAT  case 
(Figure  21)  reveal  internal  hotspots  at  both  fillet  locations 
corresponding  approximately  to  the  locations  where  the  high 
current  densities  occur.  Also  readily  observed  in  these  plots  is 
the  hot  core  along  the  inside  center  of  the  heater  element  leg 
which  is  undoubtedly  due  to  the  lack  of  radiation  cooling  on  the 
I.D.  of  the  heater  leg.  An  END-HEAT  case  was  run  at  913  kW  (6680 
A)  which  resulted  in  an  increase  of  peak  temperature  to  5622°F. 

No  experimental  data  were  available  at  this  time  to  correlate 
with  these  temperature  results. 

Figure  22  shows  the  temperature  results  for  the  K  heater 
which  was  run  only  at  the  START  HEAT  condition.  The  graphite 
base  model  results  are  shown  in  Figure  23.  The  base  results 
shown  in  Figures  23  and  24  are  for  a  slightly  higher  heater  power 
of  653  kW. 


Thermostructural  Analysis  Results — Step  3 

Although  the  heater  element  is  clamped  to  the  graphite  base 
with  set  screws  (as  shown  in  Figure  25)  both  clamped  and  free 
heater  element  leg  end-fixity  conditions  were  analyzed  to  study 
the  effects  of  end-  fixity  on  stresses.  ABAQUS  input  data  for 
the  clamped-base  J  model  are  given  in  Appendix  B.  Figure  26  and 
27  give  highly  exaggerated  deformation  plots  showing  both  the 
clamped-  and  free-end  results.  The  fact  that  the  heater  element 
leg  deflects  inward  at  the  free  end  is  probably  due  to  the 
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relatively  cool  edges  of  the  heater  element  leg,  which  do  not 
grow  by  thermal  expansion  as  much  as  the  hotter  center  portion, 
and  tend  to  pull  or  warp  the  center  toward  the  edges. 

Since  it  is  known  in  advance  that  the  heater  elements  fail 
by  brittle  fracture  at  the  top  and  bottom  fillets  and  that  the 
fracture  plane  is  approximately  perpendicular  to  the  heater 
element  longitudinal  axis,  one  might  expect  to  find  high  axial 
tensile  stresses  at  these  fillet  locations  if,  indeed,  the 
thermal  heating  stresses  are  the  cause  of  heater  failures.  The 
plots  of  axial  stress  in  Figure  28  in  fact  show  peak  tensile 
stresses  of  6896  and  7046  psi  occurring  in  the  J  heater  at  the 
top  and  bottom  fillets,  respectively,  for  START  HEAT.  The  axial 
stress  distributions  in  the  K  heater  and  the  heater  base  are 
shown  in  Figures  29  and  24.  An  assessment  of  heater  element 
failures,  based  on  these  stress  results,  is  given  in  the 
following  two  sections. 


Failure  Analysis — Stress  Based.  For  brittle-type  fractures 
in  an  isotropic  material  the  so-called  "maximum  normal  stress" 
failure  theory  can  be  used  as  the  failure  criterion.5-  This  theory 
states  that  tensile  fracture  will  occur  when  the  maximum  tensile 
principal  stress  exceeds  the  tensile  ultimate  strength  of  the 
material .  The  (algebraically)  maximum  principal  stress  plots  for 
the  J  and  K  heater  and  the  heater  base  are  shown  in  Figures  30- 
34  (Figure  35  is  discussed  on  page  16) .  Peak  tensile  principal 
stresses  at  START  HEAT  of  7000  psi*  occur  at  the  fillets  in  the  J 
heater  in  Figure  30,  but  the  K  heater  in  Figure  33  shows  a 
relatively  low  stress  at  the  fillet  and  a  relatively  low  peak 
stress  at  the  center  of  the  leg  of  4124  psi.  Table  4  lists  the 
state  of  thermal  stress  and  strain  at  each  gauss  point  in  the  top 
and  bottom  fillets  where  the  highest  principal  stress  occurs.  As 
shown  in  Table  4  by  the  "axial  tilt"  angles  of  6-19  degrees, 
these  principal  stresses  act  approximately  in  the  axial  direction 
as  expected. 


High  axial  principal  stresses  also  occur  at  the  threaded 
base  of  the  heater  element  where  it  necks  down  to  the  thinner 
wall,  and  they  are  especially  noticeable  for  the  END  HEAT  cases 
(both  clamped  and  free)  in  Figures  31  and  32.  Failures  in  this 
base  region  are  fairly  rare,  however,  indicating  that  perhaps 
there  is  some  give  at  the  base/heater  interface  which  limits 
these  stresses. 


The  principal  stress  profile  in  the  graphite  near  the 
surface  of  the  fillets  can  be  seen  in  Figures  36  and  37.  These 
figures  show  profiles  of  temperature  and  the  total  principal 

*The  axial  fillet  stresses  in  Figure  28  are  "nodal  averaged" 
results  whereas  the  principal  stresses  in  Figure  30  were  manual 
extrapolations  to  the  fillet  surfaces  of  more  accurate  gauss 
point  data.  The  use  of  two  different  extrapolation  methods 
accounts  for  the  predicted  bottom  fillet  axial  stress  being 
higher  than  the  principal  stress. 
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stress  (the  total  principal  stress  is  equal  to  the  thermal  stress 
minus  the  750  psi  compressive  stress  due  to  the  hydrostatic 
pressure  of  the  nitrogen  gas  which  surrounds  the  heater) .  The 
stress  reaches  the  maximum  level  at  the  fillet  surface  where  both 
the  temperature  and  the  strength  of  the  graphite  are  lowest.  The 
strength  of  graphite,  unlike  most  common  engineering  materials, 
actually  increases  substantially  with  increasing  temperature  up 
to  about  5000°F.  In  both  Figure  36  and  37,  the  stress  at  or  near 
the  fillet  surface  exceeds  the  indicated  4100  ksi  tensile 
strength  of  the  graphite  indicating  that  a  fracture  is  likely  to 
occur  in  the  fillets.  This  result  correlates  very  well  with  the 
failures  observed  in  the  Tunnel-9  heater  elements. 

Failure  Analysis — Strain  Based.  Table  4  lists  the  strains 
at  the  critical  locations  in  the  top  and  bottom  fillets.  Both 
the  total  strains  and  elastic  strains  are  shown.  The  elastic 
strain  (or  "strain  due  to  stress")  is  equal  to  the  total  strain 
minus  the  thermal  strain  component,  aaT,  at  the  point  in 
question.  This  elastic  strain  value  can  be  compared  to 
conventional  strain-to-failure  data  to  assess  probability  of 
failure.  Test  data  for  strain-to-failure,  taken  from  actual 
heater  element  graphite  material,  is  shown  in  Appendix  D.  The 
data  are  for  axial  and  circumferential  directions  at  an  elevated 
temperature  of  2000°F,  and  they  are  plotted  as  strain  versus 
probability-of- failure  (POF) .  A  POF  of  60  percent,  for  example, 
means  that  60  percent  of  the  test  samples  failed  at  the  indicated 
strain  level.  Profiles  of  total  principal  elastic  strain  in  the 
top  and  bottom  fillets  are  shown  in  Figures  38  and  39.  As  in  the 
stress  case,  the  maximum  strain  occurs  at  the  fillet  surface. 

The  estimated  maximum  principal  strains  indicated  in  Figures  38 
and  39  were  obtained  by  increasing  the  axial  strain  by  the 
principal  to  axial  stress  ratio  at  the  point  in  question  and  then 
subtracting  the  hydrostatic  strain*  as  shown  here: 

principal  strain  =  axial  strain /principal  stress]-  hydro,  strain 

\  axial  stress  / 

*  .00319  (5815/5471)  -  .000326 

■  .00306  in/ in  (@  Top  Fillet,  GAUSS  POINT  7) 


The  maximum  principal  elastic  strains  at  START  HEAT  are 
.00375  and  .00388  in/ in  in  the  top  and  bottom  fillets, 
respectively,  as  shown  in  Figures  38  and  39.  These  strain  levels 
correspond  to  a  POF  of  65  percent  and  75  percent,  respectively, 
which  says  65  to  75  percent  of  the  heater  elements  should  fail 
relatively  quickly,  and  that  the  remainder  should  have  relatively 
long  service  lives.  The  actual  heater  element  service  lives 
shown  earlier  in  Table  1  can  be  plotted  statistically  as  service 

*  The  hydrostatic  strain,  E^,  due  to  the  750  psi  gas  pressure  is 
computed  using:  E^  «  (l-2v)P/E  where  E^hydrostatic 
strain,  v-poisson  ratio,  P=pressure,  E«  Young's  modulus. 
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versus  percentage  of  heater  elements  which  will  fail  by  this  life 
(cumulative  frequency)  as  Figure  40  shows.  Interestingly,  the 
data  show  that  up  to  about  70  percent  of  the  heater  elements  had 
average  (20  cycles)  or  shorter  service  lives,  and  about  20 
percent  had  relatively  long  service  lives  of  50  or  more  cycles. 
This  result  agrees  well  with  the  strain  based  failure  rate 
prediction,  and  this  is  somewhat  fortuitous  considering  the  small 
statistical  population  used  and  other  factors.  This  result  shows 
that,  for  a  linear  stress  analysis  such  as  this  one,  the  strain 
based  failure  criterion  yields  more  realistic  results  than  the 
stress  based  one  which,  in  all  cases,  predicts  a  100  percent 
failure  rate  based  on  the  data  in  Appendix  D.  This  is  due  to  the 
fact  that  graphite  material  exhibits  a  strain-softening  stress- 
strain  relation  which,  if  ignored,  will  cause  the  thermal 
stresses  to  be  over-predicted. 

Table  5  summarizes  the  stress  and  strain  failure  analysis 
results.  Note  that  even  the  strain  failure  rate  is  100  percent 
for  the  high  power  (913  kW)  case  in  Table  5  owing  to  the 
substantial  increase  in  fillet  stress  with  increased  heater 
power.  It  has  been  found  that,  in  fact,  heater  element  failures 
occur  more  frequently  when  the  heater  elements  are  run  at  higher 
power  levels . 


GRAPHITE  HEATER  ELEMENT  FILLET  STRESSES 

During  the  course  of  this  investigation  various  studies  were 
performed  to: 

1.  Explain  why  high  axial  tensile  stresses  occur  at  the  top 
and  bottom  fillets. 

2.  Investigate  methods  to  reduce  high  stresses  in  the 
graphite  heater  elements,  thereby  providing  a  means  to 
extend  the  useful  service  life  of  a  heater  element. 


FUJlSt  jS£f_es?,  Breakdown 

Parametric  studies  were  performed  to  breakdown  the  bottom 
fillet  stress  into  component  parts.  Figure  41  presents  the 
results  of  these  studies.  Five  main  effects  are  identified  here 
in  order  of  importance: 

1.  Differential  Radiation  Effect  (42  percent) 

2.  Edge  Cooling  Effect  (31  percent) 

3.  Thickness  Cooling  Effect  (17  percent) 

4.  Beam  Bending  Effect  (6  percent) 

5.  Fillet  Hotspot  Effect  (4  percent) 
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Each  effect  is  described  below  briefly. 

Differential  Radiation  Effect.  The  O.D.  of  the  heater  tube 
radiates  freely  to  the  heater  vessel  liner  whereas  the  I.D.  is 
more  or  less  enclosed;  hence,  there  is  a  decreasing  temperature 
gradient  from  the  I.D.  to  the  O.D.,  tending  to  induce  axial 
tension  at  the  O.D.  and  axial  compression  at  the  I.D. 

Edge  Cooling  Effect.  The  edges  of  the  heater  element  legs 
are  more  efficiently  cooled  them  the  center  portion  because  of 
exposure  to  the  cool  nitrogen  gas  on  three  surfaces  instead  of 
just  two.  Hence  the  edges  tend  to  stay  at  a  lower  temperature 
than  the  center  and,  due  to  less  thermal  expansion,  are 
"stretched"  axially  by  the  hotter  center  and  develop  axial 
tensile  stress. 

Thickness  Cooling  Effect.  Internal  volumetric  heating  in  a 
thin  slab  with  cooling  on  the  slab  surfaces  produces  a  parabolic 
temperature  gradient.  The  core  of  the  slab  tends  to  expand  more 
than  the  outer  material  creating  a  compressive  core  and  tensile 
*  skin. 


Beam  Bending  Effect.  The  gross  thermal  expansion  of  the 
heater  element  sets  up  "beam  bending"-type  stresses  in  the  heater 
element  legs  when  the  base  is  clamped.  The  6  percent  figure 
indicated  for  this  effect  may  be  conservative  since  it  is  based 
on  a  case  where  the  heater  was  uniformly  heated  to  4000°F  which 
does  not  induce  the  more  severe  gross  warping  evident  in 
Figure  26. 

Hotspot  Effect.  Additional  thermal  expansion  of  material  at 
the  internal  hot-spots  located  near  the  fillets  tends  to  put 
additional  tensile  stress  in  the  fillets. 

The  first  three  effects  perhaps  explain  why  the  axial 
thermal  stresses  in  the  heater  element  leg  are  tensile  at  the 
O.D. ,  even  more  tensile  at  the  edges,  and  compressive  at  the 
center  core.  The  last  two  effects  may  explain  in  part  why  the 
tensile  edge  stresses  tend  to  peak  out  at  the  top  and  bottom 
fillets. 


stcssaL  Bsflasfcign  Jfe&hgfla 

Regarding  the  edge  cooling  effect,  it  was  noticed  in  the 
course  of  these  studies  that  the  temperature  in  a  heater  element 
leg  is  lowest  at  the  sharp  corners  on  the  edges  of  the  leg  (as 
shown  in  Figure  35) .  This  cooler  material  cannot  expand  as  much 
as  the  core  and,  hence,  the  comers  are  stretched  in  tension  as 
the  hot  core  expands.  It  was  found  that  rounding  off  the  sharp 
comers,  as  shown  in  Figure  35  and  42,  substantially  reduced 
fillet  stresses  on  the  order  of  20  to  30  percent. 
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Another  technique  which  may  be  useful  for  reducing  fillet 
and  edge  tensile  stresses  is  to  provide  a  "cooling  slot"  at  the 
center  of  each  heater  element  leg  to  help  cool  the  hot  core. 
Figure  43  shows  how  this  could  be  done  by  "slicing"  the  legs 
through  at  section  AA  and  rotating  the  halves  outward.  A 
prototype  design  for  such  a  heater  element  is  shown  in  Figure  44 
where  the  cooling  slot  extends  the  full  length  of  the  hot  portion 
of  the  heater  element  leg.  It  was  found  in  one  analysis  of  a 
cooling  slot  design  that  the  top  and  bottom  fillet  stresses 
decreased  by  25  percent  and  35  percent,  respectively,  as  compared 
to  the  nonslotted  heater  element. 


ASSUMPTIONS 

Various  assumptions  were  made  in  this  analysis.  For 
example,  the  use  of  one  identical  sink  temperature  for  both 
convection  and  radiation  heat  transfer  assumes  the  insulation 
liners  run  at  the  same  temperature  as  the  nitrogen  gas.  However, 
since  nitrogen  gas  is  essentially  transparent  to  heater 
radiation,  the  liner  receives  radiation  directly  from  the  heater 
element  and  probably  runs  hotter  than  the  gas.  (Preliminary 
calculations  indicate  that  the  liner  may  run  800°F  hotter.)  This 
effect  would  lead  to  a  hotter  element  and,  as  some  preliminary 
analyses  have  recently  shown,  somewhat  lower  thermal  fillet 
stresses  in  the  heater  element  (perhaps  20  percent  lower) .  on 
the  other  hand,  although  the  fillet  stresses  may  be  lower,  the 
true  ultimate  tensile  strength  of  a  given  volume  of  graphite  may 
be  significantly  lower  than  the  average  strength  used  here.  This 
is  due  to  the  actual  strength  of  a  brittle  material  (such  as 
graphite)  being  determined  by  the  distribution  of  flaws  in  the 
material.6  For  this  reason  the  tensile  strength  of  graphite  is 
highly  variable  to  such  a  point  that  the  true  strength  of  any  one 
sample  can  vary  significantly  from  the  average  strength  data. 

The  statistical  tensile  strength  plots  in  Appendix  D  were 
obtained  from  tensile  tests  on  actual  graphite  samples  which  were 
machined  from  failed  heater  element  material,  and  they  show  the 
scatter  typical  for  this  graphite.  The  statistical  strength 
variation  shown  in  Appendix  D  holds  only  for  the  specific  volume 
of  material  associated  with  the  tensile  specimen  gauge  section 
(Vol  *  0.066  cu  in) .  In  order  to  accurately  assess  failure  in  a 
graphite  part  such  as  a  heater  element,  a  sophisticated  failure 
assessment  must  be  carried  out  that  accounts  for  the  volume  of 
material  in  the  element  which  is  subjected  to  high  tensile 
stresses.  A  failure  criterion  with  these  features  requires  a 
large  body  of  statistical  material  data  which  is  currently  not 
available  for  the  heater  element.7 
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CHAPTER  4 
CONCLUSIONS 


This  report  documents  a  general  finite  element  method  for 
investigating  the  electrical,  thermal,  and  thermostructural 
performance  of  electrically  powered  heating  devices.  The  method 
was  applied  to  a  graphite  heater  element  utilized  in  the  NSWC 
Hypervelocity  Tunnel  No.  9  to  pre-heat  nitrogen  working  gas  to 
temperatures  of  3100°F.  Electric  current  density  and  temperature 
distributions  at  the  start  and  end  of  the  nitrogen  preheating 
period  are  presented.  Thermal  stresses  in  the  graphite  heater 
element  are  presented  which  predict  that  crack  initiation  of  the 
type  actually  exhibited  in  the  heater  elements  is  likely  to  occur 
at  the  top  or  bottom  fillet  of  the  J-type  heater  due  to  high 
tensile  thermal  stresses  which  occur  in  these  regions.  The 
likelihood  of  a  crack  starting  is  greater  near  the  beginning  of 
the  heating  period.  Also,  it  is  greater  with  lower  Reynold's 
Number  tunnel  conditions  since  the  relatively  low  hydrostatic 
pressure  of  nitrogen  gas  at  these  conditions  does  not  relieve  the 
tensile  fillet  stresses  enough  to  prevent  fracture.  The  results 
correctly  predict  that  crack  initiation  is  more  likely  at  higher 
heater  powers  due  to  an  increase  in  fillet  thermal  stresses  with 
increased  power.  The  stress  models  also  correctly  predict  that 
the  K-type  heater  will  tend  to  fail  at  the  bottom  fillet  because 
the  geometry  of  the  K  heater  at  the  top  fillet  effectively  limits 
the  axial  tensile  stress  in  the  top  fillet  to  relatively  low 
levels . 

Two  methods  for  reducing  fillet  stress  levels  are  presented. 
First,  the  technique  of  rounding  the  edges  of  the  heater  element 
legs  is  shown  to  reduce  fillet  stresses  by  perhaps  20  to  30 
percent.  Second,  by  providing  a  central  cooling  slot  in  each 
heater  element  leg,  the  fillet  stresses  were  reduced  by  25  to  35 
percent. 

The  results  presented  here  suggest  that  the  fillet  stresses 
generated  during  the  nitrogen  heating  period  are  a  primary  cause 
of  heater  element  fracture.  The  exact  mechanism  for  initial 
fracture  may  be  tensile  fracture  or  fatigue.  If  high  tensile 
fillet  stresses  during  heating  are  in  fact  causing  the  heater 
element  failures,  then  the  two  methods  presented  for  reducing 
fillet  and  edge  stresses,  if  implemented,  may  significantly 
increase  the  survivability  of  the  Tunnel-9  graphite  heater 
element. 
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GAS  DRIVER  VESSELS 


FIGURE  1.  SCHEMATIC  OF  TUNNEL  NO.  9  FACILITY 


FLOW 


FIGURE  2.  VERTICAL  HEATER  AND  DIAPHRAGM  SECTION 
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FIGURE  3.  HEATER  GAS  CONDITIONS  FOR  A  TYPICAL  MACH  14  WIND  TUNNEL  RUN 
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FIGURE  4.  COMPLETELY  FAILED  HEATER  -  TOP  FILLET  FIGURE  5.  COMPLETELY  FAILED  HEATER  -  BOTTOM 

LOCATION  FILLET  LOCATION 
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FIGURE  7.  PARTIAL  TOP  FILLET  CRACK 


FIGURE  8.  PARTIAL  BOTTOM  FILLET  CRACKS 
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ELECTRICAL  RESISTIVITY  (MICRO-OHM-CM) 
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TEMPERATURE 


FIGURE  9.  HEATER  GRAPHITE  ELECTRICAL  RESISTIVITY  VERSUS  TEMPERATURE 
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FIGURE  10.  J  HEATER  DRAWING 
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FIGURE  11.  K  HEATER  DRAWING 
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FIGURE  12  MACH  10  HEATER  BASE  DRAWING 
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FIGURE  13.  J  HEATER  FINITE  ELEMENT  MODEL 


FIGURE  14.  K  HEATER  FINITE  ELEMENT  MODEL 
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FIGURE  18.  CURRENT  DENSITY  CONTOURS  IN  THE  J  HEATER 
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FIGURE  20.  END  HEAT  TEMPERATURES  IN  THE  J  HEATER 
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FIGURE  26  DEFORMED  SHAPE  PLOTS  -  FLOATING  BASE 


DEFORMATION.  CLAMPED  BASE 
DEF.  MAGN  -  10X 
T-9  GRAPHITE  HEATER.  MACH  14 
5481  AMPS.  112  VDC 
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FIGURE  27.  DEFORMED  SHAPE  PLOTS  -  CLAMPED  BASE 
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FIGURE  28.  AXIAL  STRESS  IN  J  HEATER 
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FIGURE  29.  AXIAL  STRESS  IN  K  HEATER 


43/44 


FIGURE  34.  MAXIMUM  PRINCIPAL  STRESS  -  MACH-10  BASE 
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J  HEATER,  TOP  FILLET 
HYDROSTATIC  PRESSURE  -  750  PSI 
613  KWATTS  <§>  START  HEAT 
FINITE  ELEMENT  #263 


FIGURE  36.  TOTAL  PRINCIPAL  STRESS  AND  TEMPERATURE  PROFILES  EXTRAPOLATED  TO 
HEATER  SURFACE.  TOP  FILLET  LOCATION 
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J  Heater.  TOP  FILLET 
HYDROSTATIC  PRESSURE  =  750  PSI 


SURE  38.  TOTAL  PRINCIPAL  ELASTIC  STRAIN  PROFILE  EXTRAPOLATED  TO  HEATER  SURFACE, 
TOP  FILLET  LOCATION 


J  HEATER,  BOTTOM  FILLET 
HYDROSTATIC  PRESSURE  =  750  PSI 


FIGURE  39.  TOTAL  PRINCIPAL  ELASTIC  STRAIN  PROFILE  EXTRAPOLATED  TO  HEATER  SURFACE. 
BOTTOM  FILLET  LOCATION 
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CUMULATIVE  FREQUENCY 
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EFFECT 
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FIGURE  41.  BREAKDOWN  OF  THERMAL  STRESS  AT  BOTTOM  FILLET 


55 


NSWC  TR  88-146 


FIGURE  42.  FINITE  ELEMENT  MODEL  WITH  SLOT  EDGES  ROUNDED  TO  REDUCE  FILLET 
STRESSES 
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FIGURE  44.  CONCEPT  FOR  HEATER  ELEMENT  WITH  COOLING  SLOTS 
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TABLE  1.  HEATER  ELEMENT  SERVICE  LIFE  HISTORY  FOR  1985  AND  1986 

ELEMENT  LIFE 
(TUNNEL  RUNS) 


1985  FAILURES 

1986  FAILURES 

2 

10 

5 

54 

5 

9 

6 

25 

12 

2 

2 

18 

12 

7 

7 

2 

49 

87 

58 

6 

13 

1/2 

10 

5 

68 

26 

! 

21 

1 

▼ 

18  ELEMENTS.  294  RUNS 

t 

9  ELEMENTS,  228  RUNS 

AVERAGE  RUNS  PER  ELEMENT  =  19.3 
AVERAGE  ELEMENTS  PER  YEAR  =  13.5 


NSWC  TR  88-146 


59 


NSWC  TR  88-146 


< 

Q 

h- 

3 

a. 

Z 

QC 

Ul 

u. 

«/» 

Z 

< 

cc 

F- 

F- 

< 

Ui 

Z 

pi 


00 

< 


TEAT 

EMISSIVITY 

S3  1  8  S  1 

III!  !  1  !  1  f  l 

a 

z 

Ul 

1* 

BTU/FT2  HR  R 
117 

117 

117 

117 

1  1  1  1  II  1  1  1  1  1 

HEAT 

EMISSIVITY 

is  i  s  a  i 

r 

S  S  S  S  II  1  1  OS  a> 

H 

CC 

< 

5 1- 

QC 

X 

^  ifl  in 

tSKS 

□ 

CD 

100.0 

1000 

BOO 

1000 

250 

25  0 

25  0 

25  0 

100  0 

100  0 

100  0 

PRESCRIBED 

TEMP 

»■  '  '  *  *  i 

Jill  1  )  :  ;  ;  .  i 

ONVECTION 

EMPS 

END  HEAT 

R/C 

3100/3100 

-  /3100 

3100/3100 

3100/3100 

-/- 

1  i  i  1  l  ;  II. 

RAOIATION/C 
SINK  Ti 

START  HEAT 

R/C 

300/300 

-  /300 

300/300 

300/300 

-/  - 

OOS  00S 

005/ 

005/005 

005/ 

005/ 

005/ 

005/ 

005/- 

009/009 

005/009 

009/009 

REGION 

WETTED  O.D. 

WETTED  I  D 

WETTED  SLOT  EDGES 

RADIUS  CROWN 

THREADED  BASE  OO  b  ID 

WETTED  0.0 

WETTED  TOP  SURFACE 

WETTEO  BOTTOM  SURFACE 

WETTED  FLAT 

THREADED  HEATER  BASE: 

-  OD  SURFACE 

-  MATCHING  BASE  SURFACE 

THREADED  HEATER  BASE: 

-  BOTTOM  SURFACE 

-  MATCHING  BASE  SURFACE 

J/K  HEATER  OD 

J/K  HEATER  ID 

J/K  HEATER  SLOT  EDGE 

CONFIG. 

J/K 

J/K 

J/K 

J/K 

J 

BASE 

60 


FOR  913  KWATT  CASE  THESE  SURFACES  TREATED  SAME  AS  WETTED  0  D  b  I  D  OF  HEATER  LEGS 
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ESTIMATED  DATA  BASED  ON  A  SIMILAR  END  Ht AT  ANALYSIS 
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PROGRAM  CURRENT 

PROGRAM  TO  COMPUTE  TOTAL  CURRENT  IN  ABAQUS 
HEATER  MODEL  BY  SUMMING  REACTION  FLUXES  FROM 
ABAQUS  THERMAL  ANALYSIS  RESULTS.  USER  MUST  SUPPLY 
NODES  IN  THE  DATA  STATEMENT  WHICH  DEFINE  CROSS- 
SECTION  FOR  WHICH  TOTAL  CURRENT  IS  DESIRED. 

USER  MUST  UPDATE  NODEN  AND  IMPLIED  DO  LIMIT  BELOW. 

BEFORE  RUNNING  ASSIGN  THE  XXX. FI L  FILE 

TO  FORTRAN  UNIT  FOR003.  ONLY  NODES  FOR  WHICH  TEMPERATURES 
HAVE  BEEN  PRESCRIBED  CAN  BE  SPECIFIED. 

IMPLICIT  R£AL*8CA-H,0-Z) 

DOUBLE  PRECISION  ARRAY 

DIMENSION  ARRAYC 313) ,  JRRAY<2,513>  ,  NSET  C  500  > 

EQUIVALENCE  < ARRAY Cl)  ,  JRRAY (  1  , 1  )  > 

CALL  INITF 


N0DEN=29 
RFLUX=0 . 0 


DATA 

<NSET  CL) 

,L=1  , 

2?)  / 

1 

109, 

110, 

111  , 

112, 

113, 

114, 

2 

115, 

114, 

117, 

113, 

119, 

120  , 

3 

382 

,  383, 

384, 

335, 

4 

384 

,  387, 

388, 

389, 

398, 

391, 

5 

392 

,  393, 

394, 

395, 

394, 

397, 

C 

K=  1 

DO  WHILE  CK.LE.9999??) 

CALL  OBFILE  < 0 .ARRAY , JRCD) 

IF  CJRCD.NE.0)  GO  TO  100 
KEY=  JRRAY Cl, 2) 

IF  <  KEY . EQ . 21 4)  THEN 
M=  1 

DO  WHILE  CM.LE.NQDEN) 

IF  CNSETCM) .EQ. JRRAYC 1 ,3) )  THEN 
RFLUX=RFLUX+ARRAY<  4) 

END  IF 
M=M+  1 
END  DO 
END  IF 
K»K+1 
END  DO 

100  CURENTa  RFLUX*2.0 

WRITE  <4,2000)  CURENT 

2000  FORMAT  <  19X.24H  TOTAL  HEATER  CURRENT  =  F  l  0 . 3 , 1 X. ,  7HAMPERES  ) 
STOP 
END 
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PROGRAM  I2RFIL 

C  PROGRAM  TO  READ  FLUXES  AT  MAT / L  CALCULATION 
C  POINTS  <=GAUSS  PTS  FOR  DC3D20 )  FROM  ABAQUS 

C  GENERATED  XXX. FI L  FILE  AND  CONCERT  FLUXES  TO 

C  ELECTRICAL  RESISTANCE  (I2R)  HEATING  VALUES.  PROGRAM 

C  WRITES  HEATING  VALUES  TO  FILE  UNIT  032.  BEFORE 

C  RUNNING  ASSIGN  THE  XXX. FIL  FILE  TO  FORTRAN  UNIT 

C  FOR008.  PROGRAM  ASSUMES  FLUXES  ARE  CURRENT 

C  DENSITIES. 

IMPLICIT  REAL*S<  A-H , O-Z ) 

DOUBLE  PRECISION  ARRAY 
D I MENS I  ON  ARRAY  < 5 13),  JR RAY (2,513) 

EQUIVALENCE  ( ARRAY < 1 )  ,  JRRAY (1,1)) 

CALL  INITF 

C*****  SPECIFY  MATERIAL  ELECTRICAL  CONDUCTIVITY,  SIGMA 
C*****  FOP  BODY  FLUX  UNITS  OF  BTU/SEC-CU . IN . 

C*****  USE  SIGMA  UNITS  OF  (OHM-IN)**-! 

C 

SIGMA=2032.0 

DO  100  K  =  1  ,  99999 

CALL  DBF I LE( 0 , ARRAY , JRCD) 

IF< JRCD.NE.0)  GO  TO  1 10 
KEY* JRRAY  < 1,2) 

I F<  KEY . EQ . 1 )  GO  TO  10 
I F <  KEY . EQ . 28)  GO  TO  20 
GO  TO  100 
10  CONTINUE 

J EL* JRRAY < 1 ,3) 

JPNT = JRRAY ( 1,4) 

GO  TO  100 
20  CONTINUE 

RBFLUX=ARRAY  < 3 ) * ARRAY < 3 ) /S I GMA 
C 

C»****  CONVERSION  FOR  BODY  FLUX  UNITS  OF  B . T . U . /SEC-CU . I N . 
RBFLUX=RBFLUX*0 .000943451 

WRITE <32, 120)  JEL , JPNT , RBFLUX 
120  FORMAT <2I5,E10.3> 

100  CONTINUE 
110  CONTINUE 
ENDFILE  32 
STOP 
END 
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SUBROUT I NE  OFLUX  <  FLUX , TEMP , KSTEP , K I NC , T I ME , NOEL , NPT , COORDS , JLTYP) 
IMPLICIT  DOUBLE  PRECISION  <A-H,0-Z) 

DIMENSION  COORDSO),  RBFLUX< 1 000 , 27) 

DATA  KSET/0/ 

FLUXES  ARE  REAO  FROM  FILE  32  AND  LOAOED 
INTO  LOOK-UP  TABLE  RBFLUX  ON  FIRST  CALL 
OF  THIS  SUBROUTINE 

IF  < KSET . EQ . 0 )  THEN 
DO  20  1-1 ,30000 

R£AD( 32, 1000, END® 30 )  NNOEL , NNPT , BFLUX 
1000  FORMAT  <215, E10 .3) 

RBFLUX ( NNOEL , NNPT ) =BFLUX 
20  CONTINUE 

30  KSET*1 
END  IF 

C*****  DEFINE  FLUX 
C#**** 

FLUX=R8FLUX ( NOEL , NPT) 

RETURN 

END 
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PROGRAM  POWER 

*  TO  RUN  THIS  PROGRAM  YOU  MUST: 

*  <1>  ASSIGN  A  XXX. FI L  FILE  TO  A  FORTRAN  FOR008  FILE 

*  (2)  LINK/NOMAP/EXE:  POWER  POWER  ,VAXA*DUA1  :  C  ABACUS]  ABQ5LI  B/LI  B 

* 

*  THIS  PROGRAM  WILL  COMPUTE: 

• 

*  THE  TOTAL  POWER  GENERATED  IN  THE  MODEL,  TOTAL  VOLUME  OF  THE  MODEL 

*  AND  IT  WILL  SHOW  THE  AVERAGE  POWER  GENERATED  FROM  EACH  ELEMENT 

* 

*  THE  ELEMENTS  MAKE  UP  ONE  QUARTER  OF  THE  MODEL,  SO 

*  THE  TOTAL  VOLUME  AND  TOTAL  POWER  IS  MULTIPLIED  BY  4  IN  PROGRAM 

* 


IMPLICIT  REAL*8( A-H , Q-Z ) 

DOUBLE  PRECISION  ARRAY 
D I  MENS I  ON  ARRAY (513), JR RAY (2,513) 

* 

#  YOU  MAY  WANT  TO  CHANGE  THE  VALUES  DIMENSIONED  IF  YOU  HAVE  A 

*  LARGER  OR  SMALLER  MODEL 

DIMENSION  NODE<5000>,  ><(5000),  Y<5000),  2(5000) 

DIMENSION  NELE<  1000) ,NODl (  1000) ,N0D2(  1000)  ,N0D3(  1 000  >  ,N0D4( 1 000 
D I MENS I  ON  NODS (  1000)  ,  NODS (  1000),  NOD? < 1 0  0 6 )  ,  NODS (  1 0  0 0 ) 
DIMENSION  AVE( 1000) ,  NELA(1000),  VOL(1000) 

EQUIVALENCE  (ARRAY ( t )  ,  JRRAY( 1,1)) 

CALL  INITF 


L=0 

N=0 

DO  500  K  =  1 , 99999 

CALL  DBF I LE( 0 , ARRAY , JRCD) 
IF  (JRCD.NE.0)  GO,. TO  400 
KEY  =  JRRAYC 1,2) 


* 

*  KEY  1901  IS  USED  TO  OBTAIN  THE  COORDINATES  OF  EACH  NODE 

* 

IF(KEY.EQ. 1901 )  THEN 
N=N+1 

NODE(N)  =  JRRAY (1,3) 

X(N)  =  ARRAY (4) 

Y(N)  =  ARRAY ( 5) 

2 (N)  =  ARRAY (6) 

* 

*  KEY  1900  IS  USED  TO  OBTAIN  THE  NODE  POINTS  OF  EACH  ELEMENT 
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ELSE  IF  < KEY. EQ. 1900)  THEN 
L*L+  1 

NELE( L)  =  JRRAY  < 1 ,3) 

N001 <  L)  =  JRRAY  < 1  ,5) 

N0D2<  L)  =  JRRAY < 1  , 6 ) 

N003( L)  =  JRRAY ( 1,7) 

N0D4(  L)  =  JRRAY < 1  ,8) 

N0D5<L>  =  JRRAY < 1 ,9) 

N0D6<  L>  =  JRRAY < 1,10) 

NGD7' L )  »  JRRAY' 1,11) 

N0D8<  L)  =  JRRA7<1,12> 

ELSE 

T=2 
END  IF 

300  CONTINUE 

600  TOTVOL  =  0 

DO  900  I  =  1  ,L 
LEVEL  =  0 
DO  300  J  =  1  ,N 

I F< LEVEL . EQ . 8)  GO  TO  850 

IF-THEN  LOOP  CHECKS  AND  GETS  THE  COORDINATES  FOR  THE  8  NODE  PC IN 

FOR  EACH  ELEMENT,  THESE  COORDINATES  ARE  LATER  USED  TO  FIND  VOLUM 

IF  <NOOE< J> .EQ.NODt < I > )  THEN 
LEVEL=LEVEL+ 1 
PA  =J 
END  IF 

IF  (NODE< J) . EQ .NQ02< I ) )  THEN 
LEVEL  =  LEVEL  >  1 
PB  =J 
END  I F 

IF  <N0DE<J).EQ.N0D3<I)>  THEN 
LEVEL  =  LEVEL  +  1 
PC  =J 
END  I F 

IF  <NOOE<  J)  .  EGT.N0D4<  I  )  )  THEN 
LEVEL  =  LEVEL  ♦  1 
PD  =J 
END  IF 

IF  <NODE< J) .EQ.NOD5C I > >  THEN 
LEVEL  »  LEVEL  *  1 
PE  -J 
END  IF 

IF  <NOOE< J> .EQ.N0D6C I > )  THEN 
LEVEL  =»  LEVEL  ♦  1 
PF  -J 
END  IF 

I F  <NODE< J) . EQ .N0D7< I ) )  THEN 
LEVEL  -  LEVEL  +  1 
PG  »J 
END  IF 

IF  <NOOE<  J) . EQ .N0D8< I >  >  THEN 
LEVEL  =»  LEVEL  ♦  l 
PH  -J 
END  IF 

800  CONTINUE 
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*  THE  FOLLOWING  CALCULATION  WILL  DETERMINE  THE  MOL  OF  EACH  ELEMENT 

850  A J I =-X <  PA ) -X <  PB >  +X<  PC ) +XC  PD) -X< PE ) -X <  PF ) +X < PG ) +X <  PH > 

AJ4=-X<PA)-X<PB>-X<PC)-X<PD)+X<PE)+X<PF>+X<PG)+X<PH> 

A J7=-X <  PA ) +X ( PB  >  +X <  PC ) -X <  PD  > -X  < PE  >  +X ( PF ) +X < PG > -X  <  PH  > 

A J2=-Y <  PA ) -Y  <  PB  >  + Y  <  PC  > + Y  <  PD  > -Y  < PE  > -Y  C PF  ) +  Y ( PG ) + Y ( PH ) 

A J5=-Y <  PA ) -Y <  PB  > -Y  <  PC ) -Y < PD >  + Y  <  PE ) + Y <  PF ) + Y ( PG ) + Y <  PH > 
AJ8=-Y <  PA) +Y <  PB) +Y <  PC ) -Y <  PD) -Y <  PE>  +Y ( PF >  + Y < PG > -Y ( PH ) 
AJ3=-Z  <  PA) -2  <  PB) +2  <  PC) +2  <  PD) -2 ( PE) -2  <  PF) +Z  <  PG) +Z  <  PH) 
AJ6=-Z <  PA ) -2  < PB ) -2  < PC ) -2  < PD > +  2  < PE ) +  2  < PF ) +  2  < PG ) +  2  <  PH  > 

A J9=-2 <  PA) +2 ( PB) +2  < PC> -2  < PD) -2 < PE>  +2  <  PF>  +2  < PG) ~Z < PH) 


MOLA=AJ 1 *A J5*A J9+A J2*A J6*A J7+A J3*A J4*A J3-A J3*A J5*A J7 
M0LB=-AJ2*AJ4*AJ9-AJ1 *AJ6*AJS 
MOL(I)  =  0 .015.6250  *<MOLA+MOLB> 

TOTVOL  =  TOTMOL  +  MOL< I > 

900  CONTINUE 

TOTMOL  =  4*T0TM0L 


*  THIS  PART  OF  THE  PROGRAM  IS  USED  TO  OBTAIN  THE  AMERAGE  FLUX 

*  FOR  EACH  ELEMENT 

* 

SUM  =  0 
NUMB  =  1 
KB  -  0 

DO  2000  KA  =  1 ,99999 
I F<NUMB . LE . 27)  THEN 

READ <  32,1500, END-2 1 0 0 >  NEL , NPNT , FLUX 
1500  FORMAT <  2 1  5 , E 1 0 . 3 ) 

SUM  =  SUM  +  FLUX 
NUMB  =  NUMB  +  1 
ELSE 

KB  =  KB  +  1 

*  THE  AMERAGE  FLUX  IS  MULTIPLIED  BY  1054.5  TO  CONMERT  FROM 

*  BTU/SEC  TO  WATTS 

AME<  KB)  =<  SUM/27) *1054.5 
NELA< KB)  *  NEL 
SUM  =  0 
NUMB  =  1 
END  IF 

2000  CONTINUE 
2100  TOTPOW  ■  0 

«  THIS  PORTION  OF  THE  PROGRAM  IS  USED  TO  WRITE  OUT  ALL  THE  PERTINENT 

*  INFORMATION  TO  FILE  FOR054.DAT 

»  ALSO  THE  TOTAL  POWER  IS  SUMMED  UP  BELOW 

* 
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WRITE<  54 ,2150) 

2158  FORMAT ( 1  HI  , 8X , 33HTUNNEL-9  GRAPHITE  HEATER  ELEMENT  /> 

WRITE<  54 , 2200  > 

220  8  FORMAT  (IX,  6.X  ,  37HELECTR I  CAL  RES  I  STANCE  HEAT  GENERAT I  ON////  > 

WRITE <54, 2259) 

2250  FORMAT (IX, 7HELEMENT , 5X , 6 HL GLUME , 7X , SHALE  FLUX ,  7X ,  1 0 HALE  POWER  / ) 
WRITE<  54 , 2300 ) 

2308  FORMAT < IX, 10X, 1 1H< CUBIC  IN. > , 4X , 1 1 H < WATT /IN* 3) , 6X , ?H< WATTS) / > 

DO  2750  I  EL  =  1  ,  L 

POW  =  LOL< I  EL)  *  ALE< I  EL) 

WR I TE  <  54 , 250 8 ) NELE < I  EL ) , LOL < I  EL  > , ALE  < I  EL )  , POW 
2500  FORMAT < I  7 , 4X , 3<  E 1 0 . 3 , 5X ) ) 

TOTPOW  =  TOTPOW  +  POW 
2750  CONTINUE 

TOTPOW  =  4*T0TP0W 

WRITE < 54,2300) 

280 0  FORMAT < 1 X/// IX  , 2X ,  i  1 HTOTAL  POWER , 1 2X , 1 3HT0TAL  LOLUME  / > 

WRITE <54, 2900)  TOTPOW , TOTVOL 
290 0  FORMAT <  E 1 0 . 3 , IX , 5HWATTS , 7X , E 1 0 . 3 , 1 X , 9HCUB I C  I N . > 

3000  END  FILE  54 
STOP 
END 
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APPENDIX  B 
ABAQUS  INPUT  PILES 


i 


ABAQUS  THhRMAL/ELECTRIC  ANALOGUE  INPUT  FILE 
REV-J  HEATER 
•553  KWATTS 

•BOUNDARY 
INBASE, 1 1 , ,57.743 
L I  GA  ,  1  1  .  .  0  .  0 
•MATERIAL 

•CONDUCTIVITY,  TYPE* I  SO 
**  UNITS  ARE  < OHM-IN) **-1 
2932.0 

•RESTART, WRITE 
•STEP,  LINEAR 

•HEAT  TRANSFER,  STEADY  STATE 
•NODE  PRINT 

2 »  »  »  j  2 , 2  ,  ,  2 

•EL  FILE,  TEMPS,  HEAT  FLUX,  COORDS, LOADS 
•NODE  FILE 
2 , , , ,4,2, ,4 
•END  STEP 
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abaqOS  Thermal  analysis  input  file 

REV-J  GRAPHITE  HEATER 
END  HEAT 

♦MATERIAL 

♦CONDUCTIVITY,  TYPE=ISO  » 

♦♦  UN I TS»BTU/ I N-SEC-F 
9.745E-4,  79.0 

7.731E-4,  580.0  < 

5 . 7o3E-4 ,  1000.0 

4 . 537E-4 ,  1500.0 

3.773E-4,  2000.0 

3.356E-4,  2508.0 

3.055E-4,  3000.0 

2.870E-4,  3500.0 

2.685E-4,  4000.0 

2.523E-4,  4500.8 

2.291E-4,  5000.0 

♦RESTART, WRITE 

♦STEP,  CYCLE=1 00 

♦HEAT  TRANSFER,  STEADY  STATE,  TEMTOL-10 
♦RADIATE,  2ERO=-440.0 

♦♦  EMI SSIVITY= . 85 ,  STEFAN-BOLTZMAN  = 

♦♦  3.303E-15  BTU/SQ . IN . -SEC-F^*4 
CUTOUT , R5 , 3100.0,2. 807E-1 5 
CROWN , R4 ,3100.0,2. 807E-1  5 
ODDOUT , R4 , 31 00 . 0 , 2 . 807E-1 5 
♦FILM 

♦♦  UNITS  =  BTU/SEC-SQ. IN.-F 
CUTOUT  ,  F5  ,  31  0‘0 .8,  2 . 272E-0  4 
CROWN, F4,  3100.0,  2.272E-04 
ODDOUT, F6,  3100.0,  2 . 272E-0  4 
EVEN IN.F4,  3100.0,  2.272E-04 
♦DFLUX 
ALL , 8FNU 

♦EL  PRINT,  TEMPS,  HEAT  FLUX,  COORDS , LOADS 
♦NODE  PRINT 
2, ,,,2,2, ,2 

♦EL  FILE,  TEMPS,  HEAT  FLUX,  COORDS, LOADS 
♦NODE  FILE 
2 , , , ,2,2, ,2 
♦END  STEP 
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ABAQUS  STRESS  ANALYSIS  INPUT  FILE 

CLAMPED  BASE 

REU-J  GRAPHITE  HEATER 

♦MATERIAL 


♦ELASTIC, 

TYPE* 

ISOTROPIC 

1 . 31 E+06 , 

0.13, 

70.0 

1 . 34E+06 , 

0.13, 

9  32.0 

1 . 50  E+  0  6  , 

0.13, 

1  332 . 0 

1 . 65E+06 , 

0.13, 

2732 . 0 

1 . 73E+06 , 

0.13, 

3632 . 0 

1 . 69E+06 , 

0.13, 

4532 . 0 

1 . 62E+06 , 

0.13, 

5432.0 

♦EXPANSION 

1 ,  2  E  R  0=  7  0 . 0  ,  T'i 

♦  ♦  MEAN  COEFF'S  OF  EXPANSION.  REF  STRESS-FREE  TEMP=70 . 0  F 
1 . 75E-06 ,  70.0 

3.33E-04,  2500.0 
3 . 60E-06  ,  500  0.0 
♦BOUNDARY 
XZERO ,2 , ,0.0 
L I GA , 2 , ,0.0 
INBASE , i  ,,0.0 
I NBASE , 3  ,  ,0.0 
OUTBAS , 1  ,  ,  0 . 0 
OUTBAS , 3 ,  ,  0 . 0 
♦RESTART, WRITE 
♦STEP,  LINEAR 
♦STATIC 

♦  TEMPERATURE  ,  F I  LE=  1  5  ,  BSTEF  -  1  <  I  NC=  S  >  ,  ESTEP*  1  •:  I  NC=  1  ) 

♦NODE  PRINT 

•5  t  ( 

♦EL  FILE,  TEMPS,  COORDS 

o  n 

—  ,  *- 

2 

2  2  12 

♦NODE  FILE 

2  o 

*•  i  »  >»><->  *• 

♦END  STEP 
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APPENDIX  C 

FREE  CONVECTION  HEAT  TRANSFER  FILM  COEFFICIENT 
CALCULATION  FOR  MACH-14  END  HEAT  CONDITIONS 


Heater  Gas  Conditions:  Gas  =  Nitrogen 

Pressure  =  22,000  psi=  1496  atm 
Temp.  =  3100°F  (=3560°R  =  1977°K) 

First  compute  Grashof  Number  to  determine  whether  laminar  or 
turbulent  conditions  apply: 

Gr  -  g(Ts-Tg)6L3(p/u)2 


where : 


Ts=  Heater  Surface  Temperature*  5000°F  (=5460°R  =3033°K) 

Tg=»  Nitrogen  Gas  Temperature 

B  =  Coefficient  of  Expansion  =  1/Tg  =  1/3560°R 
g  *  gravity  constant  *  32  ft/ sec/sec 
L  =*  vertical  distance  from  Heater  Element  Base 
p  =  Nitrogen  Density  @  Tave 
u  =  Viscosity  of  Nitrogen  @  Tave 

TaVe  =  Average  Film  Temperature  =  (Tg-rTs)/2=  2505 °X  =  4510 °R 


Compute  density 


lave 


and  viscosity  of  nitrogen  gas: 

11  lb/cu  ft  (at  22,000  psi  from  nitrogen  table) 


u 


•ave 


6.111  x  10"5  lbm/sec-ft 


Compute 
Gr=  (32 


Grashof  Number: 

ft/sec2) ( 54 60-3 560) ( 4  ft) 3  /  11  lb/cu  ft _ 

3560°R  \  (6.111  x  10"3  lbm/sec-ft) 


Gr=*  3.686  E13 

The  following  criteria  apply: 

Laminar  Convection:  Gr  <  1.0  E+09 

Turbulent  Convection:  Gr  >  1.0  E+09 
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It  turns  out  that  the  flow  is  laminar  for  only  an  inch  or  so  from 
the  heater  element  base,  so  turbulent  flow  is  assumed  to  exist 
over  the  full  length  of  the  heater.  The  average  film  coefficient, 
h,  can  be  gotten  then  from  Figure  7-4,  Reference  C-l,  for 
vertical  cylinders. 

h  =  .021  (k/L)  (Gr*Pr)  0  • 4 

Where : 

h  =  Average  Film  Coefficient 

k  =  Thermal  Conductivity  of  gas  at  Tave 

L  =  Vertical  height  of  Heater=  4  ft 

Gr=  Grashof  Number 

Pr=  Prandtl  Number  =  Cp*u/k 

cp=  Specific  Heat  of  gas  at  Tave 

Compute  conductivity  and  specific  heat  at  average  film 
temperature : 

k  (@4510  R)  =  0.0975  BTU/hr-ft-R 
cp  "  *  0.305  BTU/lbm  R 

Then 

Pr  *  (.305) (.22)/(0. 0975)  =  0.689 
Evaluate  Film  Coefficient: 

h  =  (.021)  (0.0975)  (.689  *  3.686  E13)0*4  /4 
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APPENDIX  D 

TENSILE  STENGTH  AND  STRAIN-TO-FAILURE 
EXPERIMENTAL  DATA  FOR  HEATER  GRAPHITE 


The  Strength  Data  presented  here  was  obtained  from 
Reference  D-l  and  D-2.  The  data  are  for  2020-type  graphite 
manufactured  by  the  Stackpole  Corporation  and  YU60ST  graphite 
made  by  Ultracarbon  Corporation. 
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PROBABILITY  OF  FAILURE  (PERCENT) 


APPENDIX  D.  TENSILE  STRENGTH  AT  2000°F 
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PROBABILITY  OF  FAILURE  (PERCENT) 
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